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Abstract: 

We provide a complete set of results for the scalar 4-point function appearing in one- 
loop calculations in QCD, QED, the electroweak Standard Model and popular extensions 
thereof. Complex internal masses, which are needed for calculations involving unstable 
particles, are supported throughout, whereas complex momenta are not supported. In 
particular, for the most general, regular case we present two independent results in terms 
of 72 and 32 dilogarithms. In addition we list explicit results for all soft- and/or coUinear- 
singular cases in dimensional regularization, mass regularization, and in regularizations of 
mixed type. The exceptional case with a vanishing modified Cayley determinant, which 
hardly appears in applications, is not considered. 
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1 Introduction 

Experimental tests of tlie Standard Model of particle physics or its extensions require 
precise theoretical predictions. These are to a large extent obtained in perturbation 
theory. For an adequate investigation of most processes at the LHC or future particle 
colliders typically at least the next-to-leading order (NLO) in QCD must be known, but 
for many processes also electroweak one-loop corrections are required. Various interesting 
processes involve three or more particles in the final state and internal resonances. 

In recent years the particle theory community made a big effort to calculate processes 
for the LHC at NLO. In this context, the traditional approach of evaluating Feynman 
diagrams has been further developed (see, e.g., Ref. [[HElElll] and references therein), 
and new methods, which abandon the use of individual diagrams, have been proposed 
and applied [ El El El El E] • Details on these developments and more references to methods 
and results can, e.g., be found in Ref. [[ID]. Apart from purely numerical approaches 
[ [m [121 [13 El [15], independent of the actual method, finally at the one-loop order 
all results are expressed as linear combinations of a set of master integrals, typically 
the standard scalar one- loop integralso In four space-time dimensions [[171 [18], and in 
dimensional regularization near four dimensions [US], only the scalar 1-point, 2-point, 
3-point, and 4-point functions are needed, since higher functions can be reduced to those 
(see references in Ref. [[ID] for different variants of these decompositions). For the cal- 
culation of electroweak corrections to multiparticle processes these are required for the 
general case, i.e., with arbitrary internal masses and external momenta. For processes 
involving unstable particles these functions are also needed for complex parameters. To 
this end, we support complex internal masses while keeping momentum variables real, 
i.e. we consider the generalisation required by the so-called "complex-mass scheme" [[20], 
which is frequently employed to describe one-loop corrections to processes with particle 
resonances. 

In Ref. [118j pioneering work in the calculation of scalar one-loop integrals has been 
carried out by 't Hooft and Veltman. For the 1-point, 2-point, and 3-point functions 
compact explicit expressions are provided that are valid for complex internal masses and 
arbitrary physical momenta. For the scalar 4-point function with real masses a result in 
terms of 24 dilogarithms has been derived in Ref. [[IS], but for complex masses only a 
sketch for a result in terms of 108 dilogarithms has been given there. This result has been 
implemented into a FORTRAN code recently [[21]. In Ref. [[22] a result for the general 
scalar 4-point function has been derived involving 16 dilogarithms, but this compact result 
is only valid for real masses. 

Besides the regular scalar functions, in most calculations involving massless or nearly 
massless particles special soft- or collinear-singular cases appear. This is, in particular, the 
case for NLO calculations in QCD, where usually dimensional regularization is applied, 
or in QED and the electroweak Standard Model, where the masses of the light leptons 
are often taken to be infinitesimally small. For the 2-point function, the singular limits 
are easy to extract, and for the 3-point function the singular cases have been summarised 
in Ref. [ [23] . However, for the 4-point functions a complete set exists only in dimensional 



^ Other methods that reduce to a different set of basic integrals also rely in general on numerical 
evaluation of the basic integrals [ US] • 



regularization so far. The scalar 4-point function with vanishing internal masses and 
no or one non-zero external mass (in addition to the two non-vanishing Mandelstam 
variables) is provided in Ref . [ [251 ES] • All cases with vanishing internal masses are given 
in Refs. [ |T9l ^6\ [27] . Finally, using also various other special cases spread in the literature, 
all singular cases relevant in dimensional regularization have been collected in Ref. [[28]. 

The purpose of this paper is to provide a complete set of formulas for the scalar 
4-point function valid for complex masses and in all regular and soft- and/or collinear- 
singular cases both in dimensional and in mass regularization. In detail, we provide two 
independent results for the scalar 4-point function with arbitrary complex squared masses 
(with negative imaginary parts), which are both valid in all kinematical regions relevant 
for physical processes. One of these results is obtained similar to the method of 't Hooft 
and Veltman [ [18] and involves 72 dilogarithms, the other result follows the strategy of 
Ref. [ [22] , leading to only 32 dilogarithms. For several specific choices of complex masses 
and kinematical invariants we have numerically compared our results to the ones obtained 
with the code DOC [[21], which is based on 108 dilogarithms, and have found agreement. 
Moreover, we provide complete analytic continuations of the results of Ref. [ [22] for the 
cases involving one or more zero masses. Additionally we summarise results for the scalar 
4-point function for soft- and coUinear-singular cases in dimensional regularization, in 
mass regularization, and in regularizations of mixed type. As only restriction we require 
that particles with infinitesimally small masses do not change their mass when emitting 
a massless particle. This is the case for all applications in QCD, QED, the electroweak 
Standard Model and popular extensions thereof. 

We have derived the singular cases in two different ways. The first method is based on 
the approach of Ref. [I23]- Here, a singular integral is explicitly calculated in a suitable 
regularization scheme and then translated to all other regularization schemes by adding 
and subtracting appropriate 3-point functions. As an independent check, many integrals 
have been directly evaluated in addition, others have been taken from the literature. 
Specifically, those that do not require dimensional regularization have been derived as 
limiting cases using the methods and results of Ref. [[22], those that require dimensional 
regularization have been taken from Ref. [ [21] or from other papers, or derived as special 
cases therefrom. We have taken care to provide concise results that are valid in all relevant 
regions of phase space. 

However, we do not provide results for exceptional phase-space points, where the 
denominator in front of the dilogarithms and logarithms vanishes. These cases do only 
rarely appear in practical calculations, and can often be avoided by introducing a finite 
width for unstable particles (for examples see Refs. [[23 [30]). Either they correspond 
to Landau singularities, where the scalar 4-point function is singular, or to exceptional 
non-singular cases that practically do not show up in calculations. If needed, they can 
be derived as limiting cases of our results. We also do not provide results for complex 
momentum variables that are required for unstable external particles [[31]. 

The paper is organised as follows: In Section [21 we summarise our conventions and 
notation. Results for the general non-singular case with complex masses and for the 
corresponding massless cases are derived in Section [31 Results for all soft- and/or collinear- 
singular cases are listed in Section [H and Section [51 contains a short summary. 



Conventions and notation 

The one-loop scalar 4-point integral is defined in D space-time dimensions as 
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with po = 0, pjj = Pi — Pj, and yU is the mass scale of dimensional regularization. The 
internal squared masses mf are allowed to be complex with phases between and — 7r/2. 
If not stated otherwise, real squared masses are assumed to have an infinitesimally small 
negative imaginary part. If rrii appears in the following expressions it is defined as the 
square root of mj with positive real part. The external momenta are real, but over-lined 
squared momenta are understood to receive an infinitesimally small positive imaginary 
part, i.e., pf- = pf- + iO. Note that for D ^ 4 our normalisation deviates from the one 
usedinRefs. [[Tll|28]. 

After introducing Feynman parameters and performing the momentum integration 
(12. ip can be written as 
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where D = 4 — 2e and the integral extends over the 3-dimensional unit simplex. The 
polynomial P, which is homogeneous of degree 2 in the four variables Xq, . . . , Xs, is given 
by 
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with the modified Cayley matrix Y defined by 
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The representation (12. 2p shows the symmetry of the 4-point function resulting from per- 
mutations of the four internal lines. This symmetry is directly realized in permutations 
of the indices {0, 1, 2, 3} of the squared masses and momenta. 
In the results, the roots of quadratic equations 
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appear which we denote as — l/rj^ i and — l/rjj 2- It is easy to show that for real momenta 
and squared masses with phases between and — 7r/2 the variables Tj^ ^ vary only on the 
first Riemann sheet, i.e., they cannot become real and negative. For real masses this is 
ensured by the infinitesimal imaginary parts. The sign of the infinitesimal imaginary part 
of Vij^k resulting from the infinitesimal imaginary parts of the masses can be determined 
as 

^i^ij,k — for R-e(rjj fc) < and Im(rjj fc) = infinitesimal. 

(2.6) 
For negative t-jj^, i.e., where the infinitesimal imaginary part matters, an infinitesimal 
imaginary part of p?- yields the same result for sgn(Imrjj /;). 

Instead of using the variables ^jj /;, we sometimes employ the shorthand 
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which is — ignoring the infinitesimal imaginary part for a moment — related to the roots of 
(12. 5 p via Tiji = {mj/mi)xij and rjj^2 = ijnj/i^i)xlj^ ■ 

The 4-point function is evaluated in terms of logarithms and dilogarithms. As usual 
we choose the principal value of the logarithm such that the cut lies along the negative 
real axis. Hence one finds for a, 6, and ah not on the negative real axis 

ln(a6) - ln(a) - ln(6) = ri{a, h) (2.8) 

with 

r]{a, b) = 2m[e{- Im a) e{- Im h) e{lm{ah)) - 0(lm a) 9{\m b) e{- Im(a6))}. (2.9) 

The principal value of the dilogarithms for {1 — z) not real and negative is given by 

Li2(z) = - / — In(l-zt). (2.10) 
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More details and various useful relations for the dilogarithm can be found in Ref. [ [32] . 
Most of our results involve only functions of the type 

InUlxA, Li2U-nx,j, (2.11) 

where the variables Xj vary only on the first Riemann sheet, i.e., — vr < arg(xj) < vr. Once 
these results have been derived in the Euclidean region and for real masses they can be 
easily analytically continued by the substitutions 
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where 9 is the Heaviside step function. In most of our results we need only the analytically 
continued dilogarithm for two variables xi, X2- This simplifies to 



£^2(3^1, 3^2) = Li2 (1 - X1X2) + r]{xi,X2) In (1 - X1X2) 
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because the term proportional to the 6'-function in (I2.13P does not contribute for ra = 2 if 
the variables Xi stay on the first Riemann sheet. 

3 Regular 4-point integrals 

3.1 4-point integral for general complex masses — method similar to 't Hooft's 
and Veltman's 

In this section we evaluate the general scalar 4-point integral Dq with complex internal 
masses rrii widely following the method of Ref. [[IB]- As stated above, we implicitly 
assume that the mass squares mf possess either a non-zero negative imaginary part or 
an infinitesimally small imaginary part that is negative. Moreover, in this section we use 
the fact the the real part of mf is non-negative (a restriction that is relaxed in the next 
section). The integral is both UV and IR finite, so that we set D = 4 from the beginning, 
and (12. 2 p simphfies to 
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with the polynomial P as defined in (12. 3p . For later convenience we define the polynomials 
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Pfe(xo, Xi, X2, X3; 15) = -(1 + /3)Po + Pi + /3P2 (Xo, Xi, X2, X3), 

which are linear and homogeneous in the Xj, and the real- valued function 
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with a and (3 being real parameters. 



Before performing the first integration, we rescale the integration parameters according 
to 

X2 = (1 - a;i)^2, xs = {1- xi)^3, (3.4) 

so that the non-trivial part of the integration boundary simphfies to the 2-dimensional 
unit simplex, 

Do = t da:i /' d6 /'"^' d^ ^^ ^IzL^ . 
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Now we perform the "Euler shift" 

6 = 1/2 + a6 (3.6) 

in order to substitute variable ^ by t/2- The parameter a is chosen to obey the condition 

= Pso + "(Pso - P32 + pIo) + o^^pIo, (3-7) 

which renders to coefficient of .^| in P vanishing. At this point it is important to realize 
that the parameter a is a real quantity for all possible momentum configurations occurring 
in any scattering or decay (sub)process. For the next intermediate steps we assume 
a > and analytically continue the final result to arbitrary real a later. For a > the 
integration over ,^2 and ,^3 can be written as 

/ d^2 / d^s = / d^a / d|/2 
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Exploiting the linearity of P in 6, the ^ dependence can be isolated, 

P((l-Xi)(l-|/2-(l + «)e3),a:i,(l-a;i)(i/2 + ae3),(l-a:i)6) = (l-:^i)e3^ + 5, (3.9) 

with 

A = A{xi,y2) = Pa((l -a;i)(l -y2),xi,{l - Xi)y2,0;aj, 

B = B{x,,y2) = P{il - x,){l - y2),x,,{l - x,)y2,0). (3.10) 

The ^-integration is carried out as follows, 
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In the last integral we combined the two contributions resulting from the boundary ^3 = 
to a single integral using J^" dy2 + /!„ d?/2 = /o^d|/2- We perform the following three 
different substitutions in the three integrals of (13. lip , respectively, 

2/2 = -a + , 2/2 = 1 , 2/2 = (3.12) 

1 — Xi 1 — Xi 1 — Xi 

and introduce the shorthands Ai and Bi for the functions of Xi and ri2 that result from 
^(3^1,2/2) and -B(xi,2/2) upon making the three substitutions (13.121) that are successively 
labelled by i = 1,2,3. The functions A^ and B^ are polynomials in xi and 772, in particular, 
the Ai are of linear and the Bi of quadratic degree. This is explicit after expressing them 
in terms of P, 
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P2 = B2{xi,r]2) = P((l + a)?72,xi,l-xi - (l + a)?72,o), 
P3 = Bs{xi,ri2) = P(l -xi -r72,xi,?72,0j. (3.13) 

The Dq integral becomes 
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In this representation it is obvious that the contributions involving Ai, Bi or A2, P2 are 
zero in the special cases a = Oora = — 1, respectively. In order to eliminate the xf terms 
in the factors containing Pj in (13.141) . we again introduce Euler shifts 

T]2 = Z + PiXi (3.15) 

with three different variables /3j which correspond to the parts containing Ai and Pj. The 
following conditions render the above-mentioned coefficients of xi zero, but at the same 
time lead to real- valued values of f3i (for physical decay and scattering processes), 

= P31 + /3i(P3i - Pio + pIo) + PipIo. 
= P21 + /32(P2i - pIi + P32) + /32V32, 
= pIo + HpI^ - pIi + P20) + /3|P2o- (3.16) 

For intermediate steps we again impose the condition /3j > and perform analytical 
continuations to arbitrary real /3j at the end. Having eliminated the x\ terms in the Pj 
parts of (I3.14p . each of the three terms in curly brackets can be decomposed via partial 
fractioning in Xi, so that only terms of the form l/(axi + h) with Xi-independent terms 
a and h result. To make this step explicit, we isolate the Xi terms in Ai and Pj, 

A,{x^,z + I5,x^) = AS°)(z) + xiAf)(z), 

Pi(xi,;2 + Axi) = Bf\z)+xiBf\z) + xlBf\z). (3.17) 



The arguments of A\ and B^ are suppressed in the following. The coefficients B^ are 
fixed by the above elimination of xf terms, 
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and drop out in the subsequent manipulations. The other coefficients read 
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The partial fractioning of the integrand in (13.141) with respect to its Xi dependence yields 
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which are quadratic polynomials in z. In order to carry out the xi-integration we manip- 
ulate to double integral analogously to (13. 8p . 
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Organising the logarithms in this way, of course, requires some care concerning their cut 
structures as functions of z. For instance, we have evaluated the integrals containing only 
Ai functions as follows, 
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justified by the fact that all Tj and A] are real valued, so that the arguments of the 
two logarithms in the differences have the same imaginary part. Similar manipulations 
are possible for the integrals containing Bi functions after observing that the imaginary 
parts of A\ and B^ are z-independent and that the signs of the imaginary parts in 
the resulting logarithms do not change with varying z (i.e. no branch cuts are crossed). 
Moreover, it can be easily checked that in (13.231) the residues of the zeros of the functions 
Ci{z) add up to zero. This fact would be useful if we proceeded with the final integration 
as suggested in Ref. [HE]; however, we go a somewhat different way and exploit this fact 
at the very end in the analytical continuation in a and (3i. 

For later convenience, we factorize the quadratic polynomials Ci{z) as follows. 
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i.e. Zi^i and Zi^2 are the solutions of the quadratic equations Ci{z) = 0. The constants a, 
bi, and q are exphcitly given by 
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1 + a)Yso - aY23 - mljPab(a, a(l + /3i)j - Pa(l + a, 0, -a, 0; a) 
X Pb(l + a, 0, -a, 0; a(l + /3i)) + (1 + /3i)Pa(l + a, 0, -a, 0; a)^] , 

02 = (1 + a) [-ploPab{a, -1 - (1 + a)/32) + Pa(l, 0, -1, 0; a) 

X Pb(l, 0, -1, 0; -1 - (1 + «)/32) + Pa(l, 0, -1, 0; afP2\ , 

^2 = -iY2o - 2ml)Pab{a, -1 - (1 + a)/32) + P„(l, 0, -1, 0; a) 

X Pf,(0, 0, 1, 0; -1 - (1 + a)f32) + Pa(0, 0, 1, 0; a)Pb{l, 0, -1, 0; -1 - (1 + a)^2) 
+ 2/32P,(l,0,-l,0;a)P,(0,0,l,0;a), 

nlPab{a, -1 - (1 + a)(32) + Pa{0, 0, 1, 0; «)Pf,(o, 0, 1, 0; -1 - (1 + a)/32 



C2 



1 + a 



+ /32Pa(0,0,l,0,a) 



3J, 



as = ploPabia, (3s) - Pa(-1, 0, 1, 0; a)Pbi-l, 0, 1, 0; A 

63 = (1^20 - 2ml)Pabia, Ps) - Pa(l, 0, 0, 0; a)Pbi~l, 0, 1, 0; /^s) 

-Pa(-l,0,l,0;«)Pb(l,0,0,0;/33), 
C3 = mlPabia, /Ja) - Pa(l, 0, 0, 0; a)Pbil, 0, 0, 0; /^a), 

and obey the non-trivial relations 

det Y = b'^ — AaiCi, i = 1,2, 3, 
so that we can define 

-bi + VdetY -bi - VdetY 



-i,l 



Zi,2 



2ai ' 2ai 

In order to simplify the integrand of (13.231) . we perform the substitutions 

/■ — ^ /- — ^ ~ ^ /■ — 



(3.26) 



(3.27) 



(3.28) 



(3.29) 



in the three different blocks of integrals /q ^' dz, /l^. dz, J^ dz, respectively, leading to a 
common range of integration over C, viz. /q d^. This transforms the Zij into 



CI 



«J 



'A' ^^•^■" 1 + A' 
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Sij ^i 



J' 



(3.30) 



and integral fl3.23p becomes 
1 n 



Dn 



dC 



X 



VdetY Jo 
(C - C(,i)-' - (C - C(,2)"'] {in Ai(C, 0) - ln[(l - OMC 0) + 5i(C, 0)] } 

(C-C2,i)"'-(C-C2,2)~']{ln^(C,0)-ln52(C,0)} 
(C-C3,i)-'-(C-C3,2)^']{lnA3(C,0)-lni?3(C,0)} 
(C - Cli)-' - (C - Cr,2)"'] {in Ai(C, 1 - - lni?i(C, 1 - 0} 

(c - ^i)-^ - (c - ^2)-^; 

<{lnA2(C,l-C)-ln[(l-CM2(C,l-C) + 52(C,l-C)]} 

(C - a,i)-^ - (C - ^2)^'] {ln^(C, 1 - - lnS3(C, 1 - 0} 

(C - Ci,i)"' - (C - Ci,2)-'] {in Ai(0, C) - ln[(l - C)Ai(0, C) + 5i(0, 0] } 

(C - C2,i)"' - (C - C2,2)^'] {in ^2(0, C) - ln[CA2(0, C) + i?2(0, O] } 



(C - C3,i)-' - (C - C3,2)-'] {in ^3(0, C) - ln53(0, 0} 



Finally, we perform the projective transformation, 

c 



X 



c 



(3.31) 



(3.32) 



with 



X, 



C ■ 



«j 



^^^■,3 



«J 



i-C 



«J 



^i.j ~r Pi 



X. 



C 






J- Zij 



«J 



^ Si,j Pi + Zi,j 



X 



5»J 



■'J J 



*J 



leading to 



Dn 



Vdet Y Jo 
+ 



dx 



+ 



X — x'2 i)^^ — {x — x'2 2)^^ /2l(x) 

^ - 4,1)"^ -{x- x'^oY^ lm{x) 



[X Xi2 



I \-l 



1 - C, 



^31 (a^) 



J- ^i,j 



(3.33) 



[x x^ 2) 



/oi(x) 



X-Xi_l) 

X — X3 i)^"*^ — (x — ^32)^^ ^21 (x) 

X - Xi,i)"^ - (x - Xi,2)"^ /30(x) 

X - X2,l)""^ - (x - X2,2)""^ ^23(3^) 

X - X3,i)"^ - (X - X3,2)"^ /02(x) 



(3.34) 
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where we made use of the generic functions 

lij{x) = ln(l + x) + ln(M^ + Mfx) - ln(m,^ + YijX + m^x^) (3.35) 



InMj^ — Inm^ + ln(l + x) + Inl 1 + j^^a; j — ln(l + rij^ix) — ln(l + rij^2x), 



which depend on the constants 

(^ = 

r2 _ /\/r2 ^2 2 



M^ = -(1 + a)ml + ami + ml, 



Ml = M^ + {1 + a)p% - apl, - pi,, 

M^ = M^ + {l + a)plo-pl2, 

Ml = M^ + (1 + a)pl, - apj,. (3.36) 

In the second relation of ( 13.350 we use the variables Vij^k as defined in ( 12.50 . The validity 
of this relation for all real positive x becomes obvious after realizing that it is trivially 
fulfilled for x = and that none of the arguments of the logarithms crosses a branch cut. 
Note that the freedom to add an infinitesimally small negative imaginary part to each 
mass m,"^ should be exploited to render the imaginary part of M^ (and thus of all Mf) 
non-zero. This point is important in the cases when some masses m,i are equal, so that 
the imaginary parts entering M^ might compensate, which should be avoided. It should 
also be realized that the contributions involving the variables x^ • or x" • drop out in the 
special cases /3j = or /3j = — 1, respectively. 

The integral representation (13.34p . together with (13.35p . is well suited to perform the 
analytical continuation from a > and /3j > to any real values of a and /3j . During this 
continuation the masses rrii are kept fixed. We first vary the squared momenta pf^, p^i, 
and pI, to reach the domain where one or more /3i < 0, while keeping p|q, p^o, and ^§2 
(and thus also a, M^, Mq , M2, and M^) fixed. In this procedure, none of the arguments 
in the logarithms of the relevant lij{x) crosses a branch cut. Moreover, it does not matter 
whether one of the zeroes Xjj-, x'^j, x'-j runs over the path of the x- integration, which is 
the positive real axes, or not, because the residues at the poles at Xjj, x^^, x'-j vanish as 
explained above. This shows that the continuation to Pi < is trivial, i.e. that (13.341) is 
correct without modification. The only complication in the continuation to a < stems 
from the fact that the imaginary part of M^ (and thus of all Mf) can change sign. This 
sign change can lead to a discontinuity of ±27ri in the InM^^ term in lij{x). However, 
these possible constant terms in lij{x) compensate each other within the sum in ( 13.341) . 
because after the x-integration they are proportional to 

ln(-x:,^.) - In(-x^) - ln(-x,,) = 0. (3.37) 

This relation holds for all real values of /3j. Again no extra terms arise if any of the Xij, 
x^j, x'-j runs over the integration axis because of the vanishing residues at these poles. 
In summary, (13.341) is valid without modification for all real values of a and /3j. 

The x-integration in (13.341) can be easily carried out. Using the auxiliary integral (B.l) 
of Ref . [ I22] , the final result reads 



^0 = -TTT^ E(-l)'^H^3iK,,) + L,,ix',^,) + L,,ix',^,) - Loi(x'(,) - L,,ixl,) 
V det Y f.^1 ^ 
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with the generic function 



L2'i{x2,k 



^02 



(a^s.fc)}, 



Lij{x) 



M 



If Mf — )■ the function Li Ax) reduces to 



-[InM^ - Inm-] In(-x) + Li2(l + x) 

+ ^^2 1-^, -X - Ci2(rijA, -x) - Ci2(rij,2, -x 



2' 



(3.38) 



(3.39) 



Lij{x) j^,j2^Q 



-[InM.^ -lnm,^]ln(-x) - -ln^(-x) + Li2(l + x) 
- Ci2{rij^i, -x) - Ci2{rij^2, -x) + const., 



(3.40) 



where the x- independent term indicated as "const." drops out in the difference 
J2'i=ii~^)^~^^ ^ij i^k) ■ The final resuh (I3.38P for the Dq function with non-vanishing com- 
plex masses rrii contains 72 dilogarithms of complex arguments. 

3.2 4-point integral for general complex masses — method of Ref. [ 122] 

Following the method of Ref. [22], we first derive a result for the 4-point function 
involving 16 dilogarithms that is valid for complex masses if one of the external momenta 
squared vanishes. Then we use an identity to express the general 4-point function in terms 
of two 4-point functions of this kind. 

3.2.1 General massive 4-point integral for one vanishing external momentum squared 

We start by deriving a result for the regular scalar 4-point function with four internal 
real masses and two space-like or light-like momenta. Upon analytic continuation in the 
masses and one of the momentum squares we generalise it to the case of complex masses 
and five arbitrary real invariants, while the sixth invariant has to fulfil certain conditions, 
which hold in particular if this invariant vanishes. We follow the methods of Ref. 

In the following we need the polynomials 



and 



Q{xo,Xi,X2,X3) 



Q{xo,Xi,X2,X3] 



1 d 



+ 



d 



^20 dxQ dx 



d 



1 d 



dxi ri3 8x3 



P{xo,Xi,X2,X3) 



P{Xo,Xi,X2,X3) 



(3.41) 



(3.42) 



where r2o and ri3 are each one of the two variables defined in the corresponding equation 
(EH). 

We start from the representation (13. ip and assume that all masses are real and that 
the quantities r2o and vi^ are real and positive, which is the case if the momenta P20 and 
P31 are space-like or light-like. We apply the projective transformation 



fxA 

X2 

Us/ 



1 



/ 



1 + X + (1 - l/rn)y + (1 - l/r2o)2: 



y 

z 



\ 



(3.43) 



\x-y/ri3J 
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with the Jacobian 

d{Xi,X2,X3) ^ 1 

d{x,y,z) [l+x + {l-l/rn)y + {l-l/r2o)z]^ ^' ' 

in order to hnearise P(l — xi — a;2 — X3,xi,X2,X3) simultaneously in two variables of 
integration. The transformation of the integration region is determined by the inverse 

f x\ fxs + Xi/riA 



y 



xi - (1 - l/r2o)a;2 - x^ 



V 



Xi 
X2 



(3.45) 



/ 



This yields 



Da 



dx dy dz — — — 

-'0 Jo [P{1- z/r2o,y,z,x-y/ri3) -i€\' 



(3.46) 



where P is linear in both y and z. Note that we have included an explicit infinitesimal 
imaginary part —ie consistent with the infinitesimal imaginary part of the masses. How- 
ever, for later use the masses are still assumed to have an infinitesimal imaginary part 
which is negligible compared to ie. Performing the first integration and subsequently 
combining the resulting terms over a common denominator results in 



Dn 



dx / dy -— ; — - 

Jo Q{l,y,0,x-y/ris) 



^P(l, y,0,x- y/ris) - ie P(0, y, rao, x - y/ri^) - ie ^ 



r2o I dx I dy- 



(3.47) 



/o JO P{l,y,0,x-y/ri3) ~ieP{0,y,r2o,x -y/ris) -ie' 

By means of a partial-fraction decomposition the second integration can be done as well, 



Do 



oo ]^ 

dx — - , 
R{x) Jo 

dx — -— 
R{x) 



xri3 



dy 



g(i,o,o,x) 



Q{0,0,r2o,x) 



P(l, y,0,x- y/ris) - ie P(0, y, rao, x - y/ri^) - ie ^ 
ln[P(0, 0, r2o, x) - ie] - ln[P(l, 0, 0, x) - ie] 



- ln[P(0, xri3,r2o, 0) - ie] + ln[P(l, xr^, 0, 0) - ie] 
where the y-independent function R{x) is defined as 



R(x) = — 
r2o 



(5(1,0, 0, x) [P(0, y, rso, x - y/ri^) - ie 
Q(0, 0, r2o, x) [P(l, y,0,x- y/ru) - ie 



with 



= ax"^ + bx + c + ied = a{x — xi){x — X2) 

a = mlri3(Y23 - ^03/^20) - "^3(^12 - ^01/^20), 

b = {mlr2o - ml /r 20) {mlr 13 - ml/ris) + Y01Y23 - YqsY^ 

c = mlr2oiYoi - ^03/^3) - m^(Fi2 - ^23/^3), 

d = Yi2- Yoi/r2o - Y23/ri3 + ^03/(^20^3)- 



12, 



(3.48) 



(3.49) 



(3.50) 
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The discriminant of the quadratic form reads 

^2 _ 4oc = det F, (3.51) 

where Y is the modified Cayley matrix defined in f l2.4p . and the roots a;i^2 are given by 

,,,.i_,_,±yd;n^)^_^. ,3.52) 

Note that the infinitesimal imaginary part of the variables Xi^2 results only form the ied 
term in (I3.49p . i.e. real mass parameters enter the calculation of a;i^2 via fl3.50p as truly 
real quantities without infinitesimal imaginary parts. 

A further partial-fraction decomposition yields an integral representation valid for real 
and positive r2o,Ti3 and real masses 



a{xi - X2)Do = dx { ) 

JO \X — Xi X — X2/ 



Xi X — X2. 

The following algebraic relations with arbitrary y and z hold for R{x) and prove useful 
for the analytic continuation 

Rix) = I [P(l - z/r2o, 0, z, x) - ie] [P(l, y,0,x- y/n^) - ie] 

- [P(l, 0, 0, x) - ie] [P{1 - z/r2o, y,z,x- y/n,) - ie] } /{yz) (3.54) 
= - I Q(l - z/r2o, 0, z, x) [P(l, y,0,x- y/ris) - ie] 

- g(l, 0, 0, x) [P(l - z/r2o, y,z,x- y/ns) - ie] | /z (3.55) 
Qil,0,0,x)[Pil,y,0,x-y/rri)-ie] 

- g(l, y, 0, X - y/ns) [P(l, 0, 0, x) - ie] | /y (3.56) 
Q(l - z/r2o, 0, z, x) (5(1, y, 0, X - y/ns) 

- d [P(l - z/r2o,y, z,x- y/ris) - ie] \ . (3.57) 



In order to perform the analytic continuation, we rewrite (I3.53P as 



- / 1 1 \ /P(0,0,l,x/r2o)-i£/rio 



a(xi — X2)Dq = dx [ ) In 



\x — Xi X — X2J \ P(l, 0, 0,x) — ie ) 

~/"i3 ^^ f_}^ l__\ j^ /^i'(0,a;ri3/r2o, 1,0) - ie/rl^ 



\x — xi X — X2J \ P(l,xri3, 0, 0) — ie / 

t^l V P(l,0,0,Xfc) -ie y J(0,oo,oo/ri3) X - Xfc 
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R Rex 
Figure 1: The contour (0, i?, -R/ris) in the complex x-plane 



where the contour of the integrals in the last line results from the one shown in Figure [T] 
in the limit i? — )■ oo. Equation fl3.58p obviously coincides with f l3.53p for real masses and 



'"is, ''^20 > 0. Now we allow for complex masses and arbitrary p\^, leading to complex ris, 
while r2o is still required to be real and positive (which entails some restrictions on p\q 
and/or the masses irt^ and m\). We note that, in particular, for p\q = we can choose 
r2o = 1, so that our result holds. We here allow also squared masses with negative real 
parts, but always with negative imaginary parts, a case for which (12. ip can be directly 
analytically continued. Equation fl3.58p yields the correct analytical continuation, since: 



No cut is crossed in the logarithms of the first two integrals of fl3.58p . since all 
polynomials P have negative imaginary parts. Note that x, x/r2Q in the first integral 
and a;ri3,xri3/r2o in the second integral are real. 

If one of the poles a;i^2 crosses the integration contour in one of the first two integrals, 
it also moves into the closed contour (0, cx), cx)/ri3). The residues add up to zero 
owing to f l3.54p where one has to set iy^z) = (xri3,r2o). For real and positive xi,2 
the explicit imaginary parts in f l3.52p and fl3.58p proportional to ie are crucial for 
this cancellation, i.e. both have to be taken into account consistently (with the same 
numerical size of e). 

Finally, we have to show that there are no further discontinuities in the last term 
of (I3.58p . apart from the required ones resulting from the possible movement of Xi 2 
inside the integration contour. If Xk is outside the contour, nothing is to show, 
because the term is zero. If Xk is within the closed contour of the integral in the 
last line of (I3.58p . it can be written as Xk = u + v/ri^ with real and positive u and 
V. Using (I3.54P for x = x^, y = v and z = r2o, we find 



P(0, 0, 1, Xfe/r2o) - ie/rfo P(0, v/r2o, 1, u/r2o) - ie/r 
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(3.59) 
P(l,0,0,Xfc)-i£ P{l,v,0,u)-ie ^ ' 

For real and positive r2o, both P(0, f /r2o, 1, u/r2o) — le and P(l, f , 0, u) — ie have a 
negative imaginary part and the ratio does never cross the cut. 

If the ratio approaches the cut of the logarithm, its imaginary part can be determined 
from (I3.59P as long as ri3 has a finite imaginary part (otherwise v is not defined). 
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If ri3 is real, the last line of fl3.58p contributes only if ri3 is negative. In this case 
the imaginary part of the ratio can be determined using fl3.55p via 



P(0,0, l,Xfe/r2o) -ie/rio _ QiO,0,r2o,Xk)_ 

) 

"lo'l -I- ( m^r^.^ — rrfi I r^r^^T, 

(3.60) 



P(l,0,0,a;fe)-i5 riog(l,0,0,a;fc) 

r2o{yi2 - ^23/^3) + {mlri^ - m^/rig)^^ 



taking into account only the infinitesimal imaginary part of Xk resulting from ie. 
Obviously, the imaginary parts of the numerator and denominator of the last ex- 
pression have the same imaginary part (from Xk) near the real axis, so that the ratio 
cannot cross the cut. 

Consequently, (I3.58P provides a correct analytical continuation. We finally reemphasize 
the treatment of infinitesimal negative imaginary parts attributed to real mass squares. 
While these infinitesimal imaginary parts are made explicit in the calculation of Xk by the 
ie terms (and whenever ie appears), the infinitesimal imaginary parts of the variables rij^k 
of (12. 5p are implicitly determined from the ones in the squared masses (or equivalently 
from p?) and are negligible compared to ie, i.e., the infinitesimal imaginary part of x^Tiz 
is solely determined by the infinitesimal imaginary part of Xk- 

Using the auxiliary integral (B.l) of Ref. [122], (I3.58P can be evaluated to 

a{xi -X2)Do{pi,P2,P3,fnl,ml,ml,ml) = 

= j:{-if^'lJ:\a2(-x,/-^)-a2(-x,,ro,, 

k=i U=iL V rso / V 

- Ci2 i-Xkri3, -^^j + Ci2 i-Xkri3, roiA 

'\ 'Vl \ P(l,0,0,x,)-k / >^ 



\m, 



(3.61) 



Note that also for negative squared masses, r^ do never cross the negative real axis as long 
as the squared masses have negative (infinitesimal) imaginary parts. The result (I3.6ip 
is in particular valid for pgo ~ O5 where r2o = 1, and arbitrary real values for the other 
invariants and complex squared masses with negative imaginary parts (and positive or 
negative real parts). Owing to the permutation symmetry of the arguments, it can be 
used if any external squared momentum vanishes, p?- = 0. 

3.2.2 Result for general massive 4-point function via propagator identity 

The result of the previous subsection can be used to derive a second result for the 
general scalar 4-point integral Dq with complex internal masses. To this end, we employ 
a propagator identity of Ref. [[18] to express the general Dq in terms of two 4-point 
functions with one vanishing external momentum squared. For the latter functions the 
result (I3.6ip is applicable. Thus, we end up with a result in terms of 32 dilogarithms for 
the general Dq function with complex masses. 
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The propagator identity (2.2) of Ref. [[18] allows to express the product of two prop- 
agators in terms of a sum of two products of two propagators, 

1-7 , 7 



(<? + Pi)^ ~ ""^1 (<? + PsY — f^z {q + Py ~ M'^ L {q + Piy — rnl (g + p^y — ml 

(3.62) 

where 7 is arbitrary and 

P = 7Pi + (1 - 7)P3, 

Fixing 7 in such a way that (recall pi = pf^ and p| = p|o) 

P' = 7'pL + liplo - pIi - pIo) + P30 = (3.64) 

and inserting the identity 03.62p into (12.11) . we can express the general scalar 4-point 
function as 

n/222222 222 2\ 

Do{PlO:P21,P32^P30^P20^P31^^0^^V^2^^3) 

= lDo{0,{p2-py,pl2,plo,plo,l^pli,ml,M^,ml,ml) 

+ (1 - -^)Do{pIo,pIi, {p2 - pf, 0,p2o> (1 - ifpli^ ^0^ ^l^ ^2^ M"^) 

= iDoiplo, 7^31' (P2 - P)^ pIo, 0, P32, ml, m^, M^ m^) 

+ {l-J)Do{plo,{l-J)y,^,{p2-p)^plo,0,pl^,mlml,M\ml), (3.65) 

where 

(P2 - pY = 7(P2i - P?o) + (1 - 7)(pL - Pso), (3.66) 

and we used the symmetry of D^ mentioned after (12. 4p to arrange the arguments such 
that the 5th momentum squared vanishes. 

As mentioned above, 7 [which corresponds to — l//3i of (13.161) ] is real for all momentum 
configurations appearing in 1 — )■ n decays or 2 — )■ n scattering processes. For such 
processes, at least one of the momenta p^i, p^o, and pio is time-like. Since the interchange 
Pso ^ Pio replaces 7 by 1 — 7 and psi f-)- pso transforms 7 into I/7, the momenta can 
always be permuted such that 

< 7 < 1. (3.67) 

Then in (I3.65P the squared masses have negative imaginary parts (while the real part of 
M^ may become negative), and the result (13.610 holds. 

To be specific we select a suitable permutation of the arguments of Dq as follows: 

• If there is no t-channel-like squared momentum among the arguments of Do, the 
centre-of-mass energy squared s must be one of the arguments, and we can choose 



pIi = 


= s and 






7 = 


1 

2pii 



pIi + pIo - Pio + V (Pii + pio - pIo? - 4piiPio J , (3.68) 

where pio can be any sum of outgoing momenta and pio = pso — Psi- The choice 
( Km fulfils (Km since pl^ > (v5|j+ v^)^- 
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For one t-channel-like squared momentum t we choose p^Q = t. Since p^^ > 0, 

Pio > and pIq < (ypli — yPio) for ^W choice of hne 1 of the diagram, 7 defined 
by mM\f fulfils (13:671) again. 

If there are two t-channel-like momentum invariants ti and t2, we choose p^i = ti 
and pIq = t2- Since pIq > 0, 



7 = 7T-J- 



2P31 

fulfils (I337|) . 



P31 + P30 - Pio - V (Psi + pio - Pio)^ - 4p3iP30 , (3-69) 



In conclusion, f l3.65p together with f l3.6ip and the appropriate permutation of momenta 
and choice of 7 according to fl3.68p or fl3.69p provides a result for the general scalar 4-point 
function with non-vanishing complex masses, valid for all scattering and decay processes. 

3.3 Regular 4-point function with one or more masses vanishing 

For vanishing masses the general results of the last two sections are not directly ap- 
plicable. Instead of extracting the corresponding limits from the final results, it is more 
convenient to use the result fl3.53p which is based on Ref . [ ,22j . If certain variables diverge 
in the following formulas (e.g. for a = 0), a well-defined result is obtained by permuting 
propagator lines of the diagram, as long as the integral is non-singular. 

We first consider the case of one vanishing mass and derive a result for rri^ = and ml, 
ml, m\ non-zero and possibly complex with negative imaginary parts. It is appropriate 
to start by considering real masses and space-like or light-like momenta P20 and P31, so 
that r2o and ri3 are real and positive. In order to make use of the result (13.531) . we choose 
the solution for r2o that satisfies 

(3.70) 





-^ 0, 

r20 m2->0 


mlr2Q — ^ Yq2 


ve obtain 






P(0,0,l,x/r2o)- 


- le/rl^ _ 


Q{0,0,l,x/r2o) 



P(0,xri3/r2o,l,0)-i£/rio -i^o g(o, xri3/r2o, 1, 0) ™i^o Q(l,a;ri3, 0, 0) 

Y02 + Y23X 



,2 



>0 Yq2 + Yi2Xri3'' 



(3.71) 



where the explicit ie terms on the very l.h.s. vanish for m^ — > so that the infinitesimal 
imaginary parts are implicitly fixed by those of the internal masses contained in the Q 
functions and Yij variables. With these replacements, (I3.53P can be written as 

a{xi - X2)Do{pi,p2,P3,ml,ml,0,ml) 

d.i^ ^^i„^ ^i-"-"--) ^ 



\x — xi X — X2J \P{l,0,0,x) — ie J 

°°/^»^^/^ ^__\^^f Q(l,a;ri3,0,0) 



\x — Xi X — X2J \P{l,xri3,0,0) — ie^ 
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+ y:(-i)^inf^ 



(5(l,0,0,Xfc) \ / dx 



/. 



dx 




(0,oo,oo/ri3) X — Xk 
Yo2 + Y23X 



\x — Xi X — X2J \m'Q + Yq^x + mlx'^ — ie J 



/o Vx — xi X — X2J yiTLQ + Yoixri^i + mKxri^y — ie J 

+ E(-l)^lnf , ^02 + y23X. \ r ^io^^ 

^ \mg + yoaa^fc + mixf. - le ) J(o,oo,oo/ri3) x - Xk 

where again the contour of the last integral goes from zero to infinity on the positive real 
axis, along an arc at infinity to the direction of l/ri3 and back to the origin (see Figured] 
for R — i- 00). The quantities Xi^2 are determined from (13.491) with 

a = m^r 13F23 - "^3^12, 

b = Yo2{mlri3 - m^/ria) + Foi^23 - ^03^12, 

c = Fo2(roi - Yos/n,) - mliYu - Y^^/n^), 

d = Yi2- Y23/ns. (3.73) 

The result (13.721) is valid for real masses and real and positive ri^. The analytic contin- 
uation proceeds as for (I3.58p . We move the squared masses into the complex plane (with 
negative imaginary part) and p|^ to arbitrary real values. Thereby ri3 becomes com- 
plex and moves on the first Riemann sheet. Equation (I3.72p yields the correct analytical 
continuation, since: 

• The arguments of the logarithms in the first two integrals do not cross the cut, 
because numerator and denominator of the fractions both have negative imaginary 
parts (which in the numerator results from the masses). 

• If one of the poles xi^2 crosses the integration contour in one of the first two integrals, 
it also moves into the closed contour (0, 00, oo/ri3). The residues add up to zero 
owing to (13.561) where one has to set y = xri^,. As stressed also above, for a;i^2 that 
are close to the positive real axis the explicit ie parts in (13.521) and (13.721) are crucial 
for this cancellation. 

• If Xk is within the closed contour of the integral in the last line, it can be written as 
Xk = u + v/ri3 with real and positive u and v. Using (13.561) for x = Xk and y = v, 
we find 



P(l, 0, 0, Xk) -ie _ P(l, V, 0, u) - ie 

11^ -\- V^l-?) -I- V^o7/ -I- V^r^lll) IP 

(3.74) 



g(i,o,o,xfc) g(i,T;,o,«) 

ttiq + ■m\v'^ + TTig-u^ + Yqiv + Fosw + Yi^^uv — ie 



Y02 + Y12V + F23« 

For squared masses with negative imaginary parts both (5(1, v, 0, u) and -P(l, f , 0, u) — 
ie have a negative imaginary part and the ratio does not cross the cut. If the ar- 
gument of the logarithm approaches the cut, its imaginary part can be determined 
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from fl3.74p as long as Imri3 7^ (otherwise v is not defined). If the imaginary part 
of ri3 vanishes, the imaginary part of the ratio can be determined via [using fl3.57p ] 



P(l,0,0,a;fc) -ie Q(l,0,0,Xfc) Yqi -Yo3/ri3 + {mjrri - ml/ris)xk 



g(l,0,0,Xfc) d Yi2-Y23/ris 

taking into account only the infinitesimal imaginary part of Xk resulting from ie 
Employing the auxiliary integral (B.l) of Ref. [[22], ( I3.72p can be evaluated to 
a{xi - X2)Do{pi,P2, P3, ml, mf, 0, mj) 



(3.75) 



k=i ^ 



Y 



23 



^^2 { -Xk, T^] - 'Ci2[ -Xk, ro3,l ) -Ci2{ -Xk, ro3,2 
J'^02 



Y 



12 



- Ci2 ( -XkTu, T^] + Ci2( -XkTis, roi,i 

^02 



+ r][ -Xk,rri 



\ Qil,0,0,Xk) J \mlJ 



Ci2( -Xkri3,roi,2 
Y, 



(3.76) 



For Yq2 = ml — p^^ = and F23 7^ this result can be used after exchanging indices and 
3 (in the sense of interchanging lines and 3 in the diagram and performing the whole 
calculation with this assignment). For F02 = ^23 = 0, F02 = ^12 = 0, or Y12 = 1^23 = the 
integral becomes soft singular; these cases are covered in the next section. 



Next we consider the case of two vanishing masses and assume ml 



m 



I = 0, but 



?71q 7^ and m^l ^ 0. The other cases with two zero-masses can be obtained by obvious 
substitutionso We choose the roots such that 



ri3 ml^O 



mlris 



m^— s>0 



Yl3, 



Yol 

'^Ol.l — -^ 2' 

mf->0 TTIq 



^3^^01,2 — ^ 77^, (3.77) 



and obtain 



a{xi - X2)Do{pi,p2,P3,ml,0,0, 



TTln 



J2{-l)''^'^lCi2i-Xk,y^] - Ci2i-Xk,ro3A - Ci2 



«2 -a;fe,^03,2 



+ Ci2 ( -Xk, TT- ) - ln(-a;fc) 



for real positive 1^13, where the Xk are determined from (I3.49P with 

a = ^131^23 ^ ^^3^12 5 b = 1^02^13 + ^01^23 ~ ^ns^l' 

c = Y01Y02 - mlYi2, d = Yi2. 



■03-' 12, 



(3.78) 



(3.79) 



By going back to the corresponding integral representation (13.720 one can convince oneself 
that (I3.78P is valid for general momenta and complex masses without modification. If 



^This includes also the case with vanishing masses on opposite lines of the diagram because of 
the symmetry w.r.t. interchange of neighbouring lines, -Do(pio,P2i!-P32,P30j-P20'?'3i,'^o,"^i' "^ii"^!) = 

^0(pf0'P20>P32,P31'P21:P30>™l''"0'™2>™3)- 
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Yqi = or Yq2 = the result can be used after interchanging indices and 3, as long as 
1^23 and Yi3 do not vanish. For Yqi = ^23 = it can be written as 

a(xi - X2)Do{pi,p2, P3, ml, 0, 0, ml) 

= E(-i)'^'{-'^^2(-xfe,ro3,i) -/:^ 

- ln(-a;fc 



i2\ -Xk,ro3^2 ) ~ 2^"^ i~^k] 



Inl^Ulnf^ 



Y, 



12 



mo 



(3.80) 



which can also be used for Y02 = and Y13 = after interchanging indices and 3. Other 
degenerate cases are soft or collinear singular (see Section H]). 



For the case of three vanishing masses, consider ml 



ms = r/io 



other cases can be obtained by obvious substitutions. Using 

0, 



^03 
?^03,1 -7^„ — 2' 



ml^O m,n 



^03,2 



rrig— )-0 



0, but ml ^ 0. The 



(3.81) 



we obtain directly from (I3.78P and (I3.79P 

a(xi - X2)Do{pi, P2,P3, ml, 0, 0, 0) 



E(-i)-^£.(--.l;)-^<--.| 



fc=l 



Y 



13 



+ £z2 ( -Xk, y- ) - ln(-a;fc) 



!^01 



'<^)MM] p-«^) 



and 



O- — ^13^235 

c = Y01Y02 - mlYi2, 



b — ^02^13 + ^01^23 ~ ^03^12; 

d = Y,2. 



(3.83) 



Degenerate cases are obtained as for (I3.78p . 

Finally for four vanishing masses, (I3.82p and (I3.83P simplify to 



a{xi - X2)Do{pi,p2,P3,0,0,0,0) 



and 



E(-i 

fc=i 



ifc+i 



jCi2i -Xk, 



Y 



23 



Y, 



02 



do 



1 Y 



01 



Xk' Yi 



13 



-ln(-Xfc) 



In 



Y 



13 



Y 



12 



In 



Yr 



02 



Yr 



03 



(3.84) 



a = Y3Yn 



13-' 23, 



^02^13 + ^01^23 ^ ^03^1 



-03-' 12, 



c = FniFr 



01-' 02, 



d = Y, 



12- 



(3.85) 



Degenerate cases are mass singular. Different results for this case have been given in 
Refs. [UniET]. 
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4 Singular 4-point integrals 

Whenever mass parameters A or Aj appear, they are understood as infinitesimal. If 
not all singularities are cured by mass parameters, dimensional regularization is applied. 
Thus, in the following all formulas are valid up to orders 0{X), 0{Xi), and 0{e). For 
convenience, we define the shorthand 



c, = r(l + e)(47r)^ 



(4.1) 



and recall our graphical notation. 



pIo 




n/2 2 2 22 2 2222\ 

= ^0(Pl0>P21'P32'P30'P20'P31'"^05"^l'"^2'"^3) 



pIi 



(4.2) 



with pij = Pi — Pj. Over-lined variables are understood to receive an infinitesimally small 
imaginary part, p]j = pj^ + iO, etc., and internal squared masses other than A^ or Xf are 
allowed to be complex with phases between and — 7r/2. We categorise the singular Dq 
functions according to their numbers of collinear and soft singularities, ncoii and nsoft, 
respectively, following the classification of Kinoshita [ IMj : 

• Each of the six squared momenta pf, gives rise to a collinear singularity if the three 



variables p~j, mf, and rrij are zero or small, i.e. of 0{Xi^). Pictorially this means that 
the squared momentum difference between two propagators of small masses is small, 
e.g. if a particle with a light-like vector couples to two light internal particles. In non- 
exceptional phase-space points of scattering reactions either {^10,^32} or {^215^30} 
or {P207P31} ci-re not small or zero, so that at most four collinear singularities can be 
present at a time. By convention we take {p20'P3i} '^'^^ small in the following. 

A soft singularity arises if a light particle is exchanged between two on-shell particles. 



Formally this situation appears if rri~ is zero and p 



m" 



and p 



mi for two 



different indices j, k not equal to i. For non-exceptional phase-space points, at most 
four soft singularities can exist at a time. 

In the following we also define the regularization-scheme-independent finite part Dq 



of a singular Dq function, as described in Ref. [ 
linear combination of scalar 3-point integrals Co{i 
propagator i to a point. 

(i) Purely soft-singular cases (ncoii = 0) 



upon subtracting an appropriate 
which result from Dq upon shrinking 



There are two cases with only one soft, but no collinear singularity. Taking ttiq as the 
"soft mass" , the common finite part of these soft-singular integrals is given by 



D 



(fin) 



D, 



Co(2) 



P20 



rrin 



(4.3) 
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Therefore, knowing Dq of ( 14. 4 p one can obtain Dq of (14. 5 P from the difference of the cor- 
responding integrals Cq{2) for rriQ = and rriQ = X, which can be found in Appendix B 
of Ref. [Ii23j- The mass-regularized integral fl4.5p was taken over from (2.9) [case (i)] of 
Ref . [ |33] . The dimensionally regularized integral can be obtained by the simple replace- 
ment 

ln(A)->|i + ln(/.) + 0(e), 

which applies if only soft (but no collinear) singularities are regulated dimensionally, and 
the results are given by 




Ph 



-|-Li2(x3i)+ Yl ^^2{x:ii 

k,l=±l 



[nsoit 



3^31 



mimsiplo - ml){l - x^^) 



21n(x3i) 



2e 



-Inf- 






J 



1^21? ^32; 



V"^2-P20 
-I 

6 ' 



+ In2(x2i) + ln^(a;32) 



(4.4) 




^soft 



2:31 



mim3{plQ-ml){l - x^i) 



21n(x3i) 



-In 



m2\ 



^2 - P20 



+ H^-4i) 



+ ln^(x2i) + ln^(x32) + Li2(x: 



31; 



+ Yl '^^2(3:31 



5-^215 -^32; 



(4.5) 



k,l=±l 



where the shorthands x^- of (12. 7p are used. Note that the precise definition of Xij, which 
implies \xij\ < 1, is important for the analytical continuation. Owing to the on-shell 
conditions the masses rrii and ms are real here, but 1712 can be complex. 
The results for 1712 = are obtained by the substitution 



-21n(x3i) In 



-^ ln(x3i) 



17121^ 



ml - P20 



+ ln^(x2i) + ln^(x32) + Yl '^^2{X31, 



■^215 -^32/1 



--ln(x3i) -In 
2 vmiTTis 



Z/2 \ 



In 



k,l=±l 



Till - pI^ 



-pIo 



In 



mj - Ph 

-pIo 



^ /=±i 



(4.6) 



with z/ = yU or z/ = A, where 



X21 rni{ml-pl^) 



X32 ms{ml-pl^)' 
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(4.7) 



The diinensionally regularized integral ( 14. 4 p corresponds to "Box 16" of Ref. [[28]. 
The following three cases contain two soft singularities. Their finite part is given by 



a 



(fin) 



Da 



Co(0) + Co(2) 



P20 



(4.8) 



i.e. the full expressions for the Dq functions can be deduced from just one out of the three 
cases upon changing the regularization by subtracting and re-adding 3-point integrals (see 
Appendix B of Ref. [|23]). The purely mass-regularized case (14. lip can be read off from 
(2.13) [case (iii)] of Ref. [|33] (see also Ref. [|36]), so that the full results read 




pIi 



(n^oft = 2) 

-2x3iln(x3i) 



+ ln 



/i 



-P20 



(4.9) 




(nsoft = 2) 

-2x3iln(x3i) 
mim^ploil - xl^] 



2e K-pio 



(4.10) 




(nsoft = 2) 

_ -2x3iln(x3i) 
mimsploil - xii) 



In 



A1A2 



(4.11) 



where X31 is again defined as in (12. 7p . Owing to the on-shell conditions the masses mi 
and m^ are real. 

The dimensionally regularized integral (14. 9 p corresponds to "Box 14" of Ref. 

(ii) Cases with one colUnear singularity (ncou = 1) 



The following three cases do not contain soft singularities and have the common finite 



part 



D 



(fin) 



^0 



jpli - pIo)Co{2) + {pj, - pI)Co{3) 
iP2o - ml){pli - ml) - (pii - ml){plo - mj) ' 



(4.12) 
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The following results are valid for complex masses m2 and m^ and m^ ^ 0. For f l4.13p 
and (14.151) results for m^ = and m2 7^ result from a permutation of the arguments, 
for (14.141) we also give a result that is valid for m^ = and m2 7^ 0: 



pIo {risoit = 0) 




(Pio - ml){pj^ - ml) - (pii - mi)(pio - ml) 

[In(xi) - ln(a;2)] 



e \mi-pio 



■^(^''^ P''^^ 2U2(^^^]-2a2ix2,x^ 



+ 2 Li2 1 9 n 



+ EE(-i)'^^2(x,,xo|, 

fc=l «=3 J 



mi - pii 



(4.13) 




P30 ('^soft = 0) 



{P2o - ml){pli - ml) - (pii - ml){plo - ml) 
Ama 



X <^21n 



"^3 - Pm 



[In(xi) - ln(z2) 



ln2f!^|-4^)+ln 
m^ - P21 



2 —2 
2 / ^^3 ~ P30 



2 4 



, 2 -2 ; -2£i2(a;2,a;i^) + V V(-l)''£i2(a;fc,x/) 



21nr- ^"^^ 



(P20 - "^2)(P3i - ml) - (pii - mDiplo - ml) [ \ml - pl^ 



[In(xi) - ln(x2)] 



2 4 



-2/:^2(a;2,a;^^)-EE(-l)''^^2( 

fc=l i=3 



^k y^l ) f ' 



(4.14) 




pIo (^soft = 0) 



1 



(Pio - "^2)(P31 - "^3) - (P2I - "^2)(P30 - ml) 

Xms 



X <^21n 



ml - pIo 



[In(xi) -ln(x2)] 



+ 2Li2 p^° ^j )-2U2^^'' ^31 



' \ 2 —2 

V"^2-P21 

+ EE(-i)'^^2(x,,xo|, 

fc=l i=3 J 



2 _2 ; 2£i2(a;2,a;i^; 



(4.15) 
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with the quantities 



Xi 



XZA 



1 —1 ' 

"^2 - P20 



2^2 



^1 - Vl\ 

1 —1 ' 

mi - pi^ 



(4.16) 



ml + ml- PI2 ± V A(pi2, mi, mj) 



2ml 



2m\ 



ml + mj- PI2 T JKph, ^l ^l) 



We remark that X3 = {m2/m-i)x32 = V'^as,! = ^32,2, X4 = (m2/m3)x32 = l/r23,2 = ?^32,i, 
with X32 and r23,fc defined in (12.71) and (12. 5p . respectively. 

The results for m.2 = rn^ = are obtained by the substitutions 



2 4 



^f](-l)*^£^2(xfc,xO + 21nf 



vm^ 



k=l 1=3 
1 



Kmj- pIq 



[In(xi) - ln(x2)] 



1 



^ 2^'''^^i)"2^'''^^'^ + 



Inf ^) + Inf ^)1 [In(xi) - ln(x2)] (4.17) 



with u = \ or u = fi. 

The dimensionally regularized integral (I4.13P corresponds to "Box 13" of Ref. [128]. 
The following four cases with no or one soft singularity have the common finite part 






Dn 



Coi2) 



(pIo - pIi)Co{S) 



P20 - ^2 {P20 - ^2) {PSI - ^l) ' 

Here the mass 777.3 is real, but 7712 can be complex. 



(4.18) 




+ ln' 



Hm^ 



ml - pI^ 



l^soft 



{p\Q-m\){pl^-ml) 



+ ln 



^2 - P2I 
^\ - P20 



2e^ e 



In 



/ims 



\ml-pli 



/im3 



21nf!^^|^)lnf^ 
\mi-vinJ \mi- 



2L'.firf?)+E^4:4^.;?- 



"^2 - P20 



i=±i 



m^ - P21 rris 



~6 



3 P31 ' 
(4.19) 




+ ln' 



Xms 



mj - pii 



(n^oft = 1) 



{plo-^DiPsi-^l) 



Ph 



^ + 21nf4^^ 
L e \mi - pio 



\mi-pLJ VAVJ Kmi-piJ 



^^ \mi-pii ms 



i=±i 



TX 



(4.20) 
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Hi ^l 

-2 ' 32 

P2I "^3 



(^soft = 0) 



(pio-mi)(pii-mi) 



In 



Ama 



\mi-pi^ 



+ ln 



TTln 



pIi 



rrio 



P20 



T 



(4.21) 




('^soft = 0) 



iPw-ml){j)li 



ml) 



+ 2Li2 



pIi 
P20 ~ ^21 



In 



Amg 



Vmi-pii 



ml 



P21 



+ Y.C^2iH^,^x^ 



i=±i 



2 -2 ' ""32 

mi - pii ms 



xJ"^'' 


-pIi 


\m2 - 


-pIo 


2 ^ 




12 [' 





(4.22) 



with X32 defined as in (12. 7p . We note that i'm2/m.^)x^2 = l/'''23,i = ''"32,2! ("^2/''^3)a^32 = 
1/^^23,2 = ^32,1, with Tij^k defined in ([23]). 

The above results are also valid for m,2 = 0, even if P32 = fn^. For m2 — )• the variable 
3:32 behaves as 



fn2 +1 _ , 2\ "^2 _i 

_ ^32 ~ ^[''^2) J _ ^32 

?7l3 ?7l3 



"^1 - Ph 

2 

mi 



+ 0{ml). 



(4.23) 



The dimensionally regularized integral (I4.19p can be found in (A. 4) of Ref. [|35] and 
corresponds to "Box 12" of Ref. [^28]. Integral (14.201) is case (ia) of Ref. [|33] (modulo 
the translation from mass to dimensional regularization which is trivial for the purely 
soft-singular case) and corresponds to (14. 4p for small mi. 

The following three cases with up to two soft singularities have the common finite part 

Co(2) C„(3) ^^^2^1 



D 



(fin) 



Dn 



P20 - ml 



P31 - rnl ' 



Here m,2 and ms have to be real quantities, 




In (X32) 



P32 



I'^soft 



1 



(P20-"^2)(P31-"^3 

+ ^lnf^^ 
e \mi-pii 



e2 e Vmi-pio 



2 In 



/im2 



/im3 



rrin 



pIo 



J \mi - 



3 P31 
(4.25) 
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ph 



(risoit = I] 



iP2o - mi){pii - mi) [ e 



'-ln( 



Xms 



Vmi-pii 



+ 2 In 



f- 



/im2 



\rni-pio 



H 



Ams 



ml - pii 



TT 



- In (X32) - — 



(4.26) 




(^soft = 0) 



1 



{p1o-ml){p: 



31 



TTln 



n 



- \n\xs2) - - 



2 In 



Xm2 



ml 



P20 



In 



Amg 



?71q 



P3I 



(4.27) 



with again X32 defined as in {\'2.7\i . 

The dimensionally regularized integral (I4.25P can be found in (A. 4) of Ref. [|35] and 
corresponds to "Box 11" of Ref. [[28]. Integral (14.261) can be calculated from case (i) 
of Ref. [|33] for m,Q = m\ and mi small (with the trivial translation to dimensional 
regularization) or from (14. 4p for p\^ = m\ and small mi. 



(in) Cases with two collinear singularities on adjacent legs {n, 



coll 



Here we distinguish six different cases with no or one soft singularity which have the 
following common finite part 



D, 



(fin) 



Du- 



(P31 - Pl2)Co{0) + {pli - pI,)Co{2) Co(3) 



P2o{pIi - ml) 



P31 



ml' 



(4.28) 



The mass m,^ can be complex, 




('^soft = 1) 



Ploip"^ 



'31 



ml) 



Q Q 



In 



+ ln 



ml - pIo 

ml - pli 



H 



-P20 



+ ln 



+ ln 



rrin 



Pm 



ml 



pIi 



lln= 

2 

2 
mi - 



/^^ 



-P20 



PI2 



ml 



2 Li. 



P32 



P31 



ml - pli 



\ml-pli 



Pm 



mq 



--In^ 
2 

pIo 

P32 -ml' ml- pIq 



m,n 



Ph 



ml 



pIi 



TT 

T 



(4.29) 
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+ ln 



(risoft = 1) 



_^ _L ^ 



Plomi - rnl) \2e'^ e 



In 



/iA 
-P20 



+ ln 






_ln2(^]+iln2 



A 



/"^ 



-P20 



/^^ 



-P20 



In 



ml - P32 
mi - pii 



+ ln 



A2 



-P20 



^^/mi-pio 



\mi-pii 



-In^ 
2 






P20 



?71q 



P32 -mf ml- pIq 



TT 

y 



rraS — 



Ph 



rrin 



pIo 



(4.30) 




A2 



+ ln 



-P20 



In 



rnl - pIo 



rrin 



P31 



l^soft 



PloiP'^ 



31 
1 



rriq 



+ ln 



+ Ci2{-2 
■-P32 



Plo 



rrin 



1 ' 2 



pIo 



-P20 



In^ 



In 






^lnfMi)+ln 



-P20 



In^ 



-P20 



"^i - pi2 



mq 



pIi 



^' 



-P20 



In 



^3 - P30 



A1A2 
-P20 



(4.31) 




VP32 -ml' mj- plo 



{risoft = 0) 



Ploip'^ 



'31 



+ ln 



r/iq 



?71q 



In 



A^ 



-P20 



pIo''' 



?71q 



pIi 



7r2~ 



ln("^^ 



rriq 



Pi2 



P31 



(4.32) 




+ jCt2[- 



pIi 
Plo 



m% 



Ph 



1 1 1 
m^ m^ 



ph 



(^soft 



0) 



P20(P^ 



'31 



ml) 



In 



+ ln 



^l - pIo 

ml - pii 



+ 



vr 



A^ 



-P20 



+ ln 



?Tlq 



ph 



ml 



P31 



V"^3-P30 



(4.33) 



30 




-In 

2 



+ Cl2 



2 —2 
2 / "^3 ~ P32 



TTlo 



P20 



pIo 



2' 



2 Li. 



m% 



Ph-^3 m 



pIo 



{risoit = 0) 






-In^ 



A^ 



-P20 



Ph' 



pIo' 



\-vinJ V \mi-vi.J \mi-vi.J 



\m-pii 



ml - pii , 



P32 P31 
ml - pli 

--]■ 

6 



\mi-pii 



(4.34) 



The dimensionally regularized integral fl4.29p was calculated in Ref. [ EZ] and corresponds 
to "Box 8" of Ref. [ |28] . Integral f l4.30p can be obtained from fl4.19p for 7713 small. Integral 
fl4.3ip is case (ib) of Ref. [ |33] (modulo the trivial translation to dimensional regulariza- 
tion) and corresponds to (14.201) for small m^. 

The following eight different cases with up to two soft singularities have the common 
finite part 
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The dimensionally regularized integral f l4.36p . e.g., appeared in the NLO QCD correction 
to pp — )• ttH [!35] [see Eq. (A. 4) there] and corresponds to "Box 7" of Ref. [128]. Integral 
fl4.37p can be obtained from fl4.25p for m2 small. Integral f l4.38p can be obtained from 
f l4.19p for P32 = m| and m^, small. Integral fl4.39p is case (iic) of Ref. [|33] (modulo 
translation to dimensional regularization) and corresponds to (14. 4 p with p\^ = ml small 
and 7712 = 0. Similar comments apply to the next two cases. Integral (I4.40p results from 
case (ic) of Ref. [|33] and corresponds to (14. 4 p with pji = and m2 = mi both small. 
Integral (I4.4ip coincides with case (i) of Ref. [ |33] for m^ = m\ and both mi and m4 small 
and corresponds to (14. 4p with p^g = ''^2 ^^d both mi and ms small. 

The following six different cases with up to three soft singularities have the common 
finite part 
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Again the mass ma has to be a real quantity here, 
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The integral (I4.45P was given in Eq. (A. 4) of Ref. [[38] and corresponds to "Box 6" of 
Ref. [[28]. Integral (I4.46P can be obtained from (14.251) for pig = "^2 ^'^d m^ small. Integral 
(I4.47P is case (iiia) of Ref. [ [33] (again to be translated to dimensional regularization, also 
in the next two cases) and corresponds to (14. 9 p for small m\. Integral (I4.48P results from 
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(iv) Cases with two colliTiear siTigularities on opposite legs (ricoii = 2) 

We distinguish seven different cases which do not contain a soft singularity and have 
the following common finite part 
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The integral (I4.52p can be found in Refs. [ [191 ISSl |39] and corresponds to "Box 3" of 
Ref. [[25. Integrals (I4.53p - p.58p can be obtained from (I4.13p - p.l5p with appropriate 
substitutions. 



(v) Cases with three collinear singularities {ricow = 3) 

There are ten different cases with up to two soft singularities and the following common 
finite part 
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The integral fl4.60p can be found in Refs. [ [T21 |2S1 El], corresponds to "Box 2" of Ref. 
and had aheady been calculated in Refs. [[211 125]. Integral f l4.6ip can be obtained from 
(A. 4) of Ref. [|35] and corresponds to (14.361) with small m^. Integral (I4.62p corresponds to 
(I4.25P with both 1112 and m^ small. Integral (I4.63P can be obtained from (I4.19P for m\ = 0, 
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3' 



and ml small. Integral (I4.64p is a particular case of (I4.19p . Integrals dH 



Ph 

and (I4.66P are cases (iid) and (id) of Ref. [|33j (translated to dimensional regularization) , 

respectively, and can also be derived from (I4.39P and (I4.40p for small ma, respectively. 

(vi) Cases with four colli7iear si7igularities (ncou = 4) 

There are six different cases with up to four soft singularities and the following common 
finite part 
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The integral (HTTj) is given in (4.25) of Ref. [[H], corresponds to "Box 1" of Ref. [[28] 
and had aheady been calculated in Refs. [[231 125]. Integral (I4.72p can be obtained from 
Ref. [[38] and corresponds to (14.45^ for small 777,3. Integral (14.73^ is a special case of (14.25^ 
for p'^2 = ^^d small masses. Integral (14.741) is case (iii b) of Ref. [ [33] (in dimensional 
regularization) and corresponds to (14.471) for small 7773. Integral (I4.75P is case (ie) of 
Ref. [[33] (in dimensional regularization) and corresponds to a special case of (14.41) . for 
vanishing pgi and P32 and if again all masses are small. 



5 Summary 

The calculation of scalar one-loop integrals is at the heart of any perturbative calcula- 
tion of radiative corrections at the one-loop level, independent of the method employed. 
In diagrammatic calculations, which do not rely on numerical integration of loop integrals, 
tensor integrals as well as any one-loop integral with more than four external particles 
(legs) are algebraically reduced to standard scalar integrals with at most four propaga- 
tors. Recently so-called unitarity-based methods became more and more popular at the 
one-loop level, i.e. methods that abandon the direct use of Feynman diagrams, but also 
those methods need the standard scalar integrals with up to four propagators as input. 
The purpose of this paper is to close some gaps in the literature and to deliver a complete 
set of 4-point ("box") scalar integrals that appear in any kind of one-loop calculation. 
For the lower-point integrals exhaustive lists of results already exist. 
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The considered scalar box integrals are naturally categorised into finite, regular cases 
and singular cases that involve soft and/or coUinear divergences. For the most general reg- 
ular case with complex internal masses, which appear in applications that involve unstable 
intermediate particles, we have presented two new analytical expressions that involve less 
(72 and 32) dilogarithmic functions than a recently published independent result of other 
authors. In addition we have listed a complete set of soft- and/or collinear-singular cases 
which appear in the calculation of strong or electroweak radiative corrections, i.e. we sup- 
port both dimensional regularization and regularizations schemes that employ small mass 
parameters. Mixed regularizations of these types are supported as well, and all results are 
valid for complex masses of internal particles. Many of those results have not yet been 
published elsewhere. 

In combination with available results on scalar 1-, 2-, and 3-point integrals now all 
relevant standard scalar integrals for one-loop calculations in QCD, QED, the electroweak 
Standard Model and typical extensions thereof are known for all not too exotic (mass or 
dimensional) regularization schemes, covering also internal unstable particles appearing 
as resonances. 

Acknowledgements 

This work is supported in part by the European Community's Marie-Curie Research 
Training Network under contract MRTN-CT-2006-035505 "Tools and Precision Calcula- 
tions for Physics Discoveries at Colliders" . 

References 



[1] 
[2] 

[3] 

[4] 
[5] 

[6] 

[7] 

[8] 
[9] 



F. del Aguila and R. Pittau, JHEP 0407 (2004) 017 |hep-ph/0404T20 



T. Binoth, J. P. Guillet, G Heinrich, E. Pilon and C. Schubert, JHEP 0510 (2005) 
015 |arXiv:hep-ph/0504267| . 



R. K. E lhs, W. T. Giele and G. Zanderighi, Phys. Rev. D 73 (2006) 014027 
arXiv:hep-ph/0508308|. 



A. Denner and S. Dittmaier, Nucl. Phys. B 734 (2006) 62 | arXiv:hep-ph/0509141 
G. Ossola, C. G. Papadopoulos and R. Pittau, Nucl. Phys. B 763, 147 (2007) 



arXiv:hep-p h/0609007| . 

Z. Bern, L. J. Dixon and D. A. Kosower, Annals Phys. 322, 1587 (2007) 
|arXiv:0704.2798l [hep-ph]]. 

R. K. Elhs, W. T. Giele and Z. Kunszt, JHEP 0803, 003 (2008) [arXiv:0708.2398l 

[hep-ph]]. 

C. F. Berger et al, Phys. Rev. D 78 (2008) 036003 [arXiv:0803.4180l [hep-ph]]. 



A. van Hameren, JHEP 0907 (2009) 088 [arXiv: 0905. 10051 [hep-ph]]. 



41 



C. Buttar et al. [QCD, EW, and Higgs Working Group], |arXiv:hep-ph/0604120 



[lo: 
[11 

[12 
[13 

[14 

[15' 
[16 

[17 
[18 
[19 

[20 

[21 

[22 
[23 
[24 
[25 
[26 
[27; 
[28 
[29 



[30] F. Boudjema and L. D. Ninh, Phys. Rev. D 78 (2008) 093005 |arXiv: 0806. 14981 [hep- 

ph]]. 

42 



Z. Bern et al. [NLO Multileg Working Group], arXiv:0803.0494 [hep-ph]; 

T. Binoth et al. [SM and NLO Multileg Working Group], larXiv: 1003. 1241. [hep-ph]. 

A. Ferrogha, M. Passera, G. Passarino and S. Uccirati, Nucl. Phys. B 650, 162 (2003) 
|hep-ph/0209219| . 

Z. Nagy and D. E. Soper, JHEP 0309 (2003) 055 |hep-ph/03 08127|. 

Y. Kurihara and T. Kaneko, Comput. Phys. Commun. 174 (2006) 530 
[ arXiv:hep-ph/0503003 | . 

S. Catani, T. Gleisberg, F. Krauss, G. Rodrigo and J. C. Winter, JHEP 0809 (2008) 
065 [ arXiv:0804.3170. [hep-ph]]. 

W. Kilian and T. Kleinschmidt, larXiv:0912.3495] [hep-ph]. 

T. Binoth, G. Heinrich and N. Kauer, Nucl. Phys. B 654 (2003) 277 
|hep-ph/0210023j . 

D. B. Melrose, Nuovo Cimento XL A (1965) 181. 

G. 't Hooft and M. J. Vehman, Nucl. Phys. B 153 (1979) 365. 

Z. Ber n, L. J. Dixon and D. A. Kosower, Nucl. Phys. B 412 (1994) 751 
|hep-ph/9306240]. 

A. Denner, S. Dit tmaier, M. Roth and L. H. Wieders, Nucl. Phys. B 724 (2005) 247 
|hep-ph/0505042] . 

T. N. Dao and D. N. Le, Comput. Phys. Commun. 180 (2009) 2258 [arXiv:0902.0325l 
[hep-ph]]. 

A. Denner, U. Nierste and R. Scharf, Nucl. Phys. B 367 (1991) 637. 



S. Dittmaier, Nucl. Phys. B 675 (2003) 447 |hep-ph/ 0308246l. 

K. Fabricius and I. Schmitt, Z. Phys. C 3 (1979) 51. 

S. Papadopoulos, A. P. Contogouris and J. Ralston, Phys. Rev. D 25 (1982) 2218. 



G. Duplancic and B. Nizic, Eur. Phys. J. C 20 (2001) 357 |arXiv:hep-ph /0006249] . 

G. Duplancic and B. Nizic, Eur. Phys. J. C 24 (2002) 385 |arXiv:hep-pii/0201306] . 

R. K. Ellis and G. Zanderighi, JHEP 0802 (2008) 002 [ arXiv:0712.185r [hep-ph]]. 

A. Denner, S. Dittm aier and T. Hahn, Phys. Rev. D 56 (1997) 117 
arXiv:hep-ph/9612390]. 



[31] G. Passarino, C. Sturm and S. Uccirati, Nucl. Phys. B 834 (2010) 77 [arXiv: 100 1.33601 
[hep-ph]]. 

[32] R. Lewin, Polylogarithms and associated functions (North-Holland, Amsterdam, 
1981). 

[33] W. Beenakker and A. Denner, Nucl. Phys. B 338 (1990) 349. 

[34] T. Kinoshita, J. Math. Phys. 3 (1962) 650. 

[35] W. Beenakker et al, Nucl. Phys. B 653 (2003) 151 [ hep-ph/02li352] . 

[36] P. Van Nieuwenhuizen, Nucl. Phys. B 28 (1971) 429; 

F. A. Berends, K. J. F. Gaemers and R. Gastmans, Nucl. Phys. B 63 (1973) 381. 

[37] E. L. Berger, M. Klasen and T. M. P. Tait, Phys. Rev. D 62, 095014 (2000) 
|arXiv:hep-ph/0005196 | . 

[38] W. Beenakker, H. Kuijf, W. L. van Neerven and J. Smith, Phys. Rev. D 40 (1989) 
54. 

[39] R. K. Ellis, D. A. Ross and A. E. Terrano, Nucl. Phys. B 178, 421 (1981). 



43 



